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Abstract 

We obtain model independent bounds for the form factors which arise 
in semileptonic B ^ tt decays. To this end we derive a theoretical 
restriction for possible combinations of the value of the form factor 
and its derivatives at the end point region. These restrictions are 
then used as an input in a — slightly modified — general formalism, 
which has already been deduced and applied in earlier publications. 
Our results can be used to set constraints on the form factors which 
have to be obeyed by all model dependent parametrizations. 



1 Introduction 



One important parameter for the determination of the unitarity triangle of 
the CKM matrix is Vub- While its phase can only be determined from non- 
leptonic decays, its magnitude will be determined from semileptonic pro- 
cesses. Various methods have been proposed to extract Vub from both in- 
clusive as well as exclusive semileptonic decays, and the final numbers after 
the era of the B factories will probably be obtained from a mixture of both 
exclusive and inclusive methods. 

The disadvantage of exclusive decays is their dependence on form fac- 
tors which are nonperturbative quantities. Modelling these form factors 
necessarily introduces some uncontrollable systematic uncertainty into the 
determination of Vub and hence a model independent method is desirable. 
At present the values of Vub extracted from exclusive semileptonic S — >■ tt 
decays using lattice QCD, QCD sum rules and quark models range from 
\Vub\ = (2.5 — 4.5) X 10^'^. Clearly this is not a satisfactory situation, in 
particular in view of the data to be expected from the B factories. Although 
there have been recent efforts to update existing models (see e.g. 
there is still room for some improvement. 

There are constraints on form factors originating from their analyticity 
properties and the unitarity of the underlying theory. The ideas how to 
implement these constraints are in fact quite old and have recently 

attracted renewed attention. In particular, combining these unitarity bounds 



with lattice data [|10[ has lead to relatively tight and model independent 
bounds on the form factors oi B ^ n transitions. 

In the present paper we improve the application of unitarity bounds in 
such a way that not only points at which the form factor is known (e.g. from 
lattice data or heavy quark relations) can be included, but also the slope and 
higher derivatives at some point, which could be known for instance from 
sum rule considerations. At the point of maximal momentum transfer the 
derivatives are correlated with the value of the form factor due to analyticity 
and we shall discuss the restrictions obtained from this in some detail. 

A similar improvement of the unitarity bounds has been proposed by 



Boyd et al. 12], where the form factors for B jriu are expanded in a 
cleverly chosen conformal variable z and some of the relations we use in the 
present paper appear in a similar form in these references. 

In the next section we shall summarize the formalism as described in . 
The improved bounds are considered in section ^ where we generalize the 
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method to include known slopes and higher derivatives. We make use of the 
correlation between the derivatives and the value of the form factor taken at 
Qmax to tighten the bounds even without using any additional physics input. 
In order to get useful bounds some physics input is necessary and we chose 
to use the chiral limit which is valid close to the end point. This has a strong 
effect since any knowledge of the form factors in the end point region tightens 
the bounds significantly. With this input we study numerical examples. 



2 Bounds on Form Factors 

For later use we summarize in this section the formalism how to derive bounds 
on the form factors using only perturbative QCD, unitarity and the analyti- 
city of the form factors in the complex plane. 

We choose to describe the hadronic matrix element of the semileptonic 
Tr~^i~Pi decays with the following two form factors: 

where = m%, — m^, — u'j^b and q — p — p' ■ In this notation 
q^ runs from = ml to q^ax = — m)^. Throughout this paper we 
will neglect the lepton masses and therefore set q^^^ — 0- Furthermore, this 
special choice of decomposition into Lorentz vectors leads to form factors 
which have to satisfy the kinematical constraint 

r (0) = r(o). (2) 

The form factors are real functions of a real variable but it is convenient to 
think of them as analytic functions in the complex g^-plane. 

To derive bounds on f^{q^) and f'^{q^) we consider the two-point function 

W = ij d^x e*«-^(0|T{\/^(a;)\/^t(0)}|0) 

^ -{g'^q' - gV)nT(g') + q'q'IlLiq'), (3) 

where = u^y^^b as above and with Urp/Liq"^) corresponding to the propaga- 
tion of a J-^ = l~/0+ particle. This two-point function can be (and has been) 
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reliable evaluated in the deep euclidean region where —q^ = Q"^ ^ ^qcd; 
at an energy scale where perturbative QCD is appliccable. 

Including a sum over all possible intermediate states F we get the following 
result for the imaginary part of (^: 

= ^E(2vrr^^'H9-P7)(o|^'^(o)|r)(r|\/^t(o)|o). (4) 

^ r 

We therefore get an equation for the spectral functions IyuHt / liq^) (i-e. for 
the absorptive parts of Hr/Liq^)) which we can relate to the real parts using 
the substracted dispersion relations 



and 




(-QW)) = -/„ (6) 



We now restrict ourselves to include only contributions of the B* and the 
B-K states in the sum over all intermediate states in (^. It is possible to 
discard all other intermediate states because @ is a sum of positive terms if 
the indices are treated symmetrically. Using isospin and crossing symmetry 
we can use (|1|) to express (||) by the two form factors. Projecting out now 
the transversal/longitudinal parts, one gets the inequalities 

and 

imnT(t) > ^(^) ^{t-^l') 

where t = and t± = (M ± m)^. 

It is now possible to get bounds on the form factors if one inserts (^,|D 
in Since the l.h.s. of (||,|D can be calculated for ^ ^qcd using 
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perturbative QCD one gets inequalities which restrict the form factors. These 
inequahties take the form (in shorthand notation) 



1 r°° 

jm>- dt k{t,Q^)\f{t)\' (9) 

where J{Q^) denotes the QCD input (i.e. the perturbative calculation of 
Xt/l) including the _B*-resonance in case of /"*" and kit^Q"^) is a known 
kinematical function. The exact value/structure of J{Q^) and k{t,Q'^) does 
of course depend on the form factor under examination. 

To translate the inequality (^ in constraints on the form factor for values 
of t in the range [0,t_] (which is the kinematical region of physical interest), 
we map the complex t-plane into the unit disc with the conformal transfor- 
mation 



l + z lU-t ^^^^ 



1-z \ t+-t_ 
so that (H) becomes 

JiQ')>[ ^mz,Q')f{z)\' (11) 

J\z\=i Zmz 

with f{z)=f{t{z)). Here we have used the fact, that k(t,Q^) is a positive 

definite quantity so that 0(z, Q^) = ^k{t{z), Q'^) times the squareroot of the 
Jacobian of the transformation. 

The value of the form factor f{z) for any point z{t) is accessible by 
defining an inner product 

and by considering the product {gt\4>f) where 

= r^i(tj^' 

so that f{z{t)) has no poles in the range [0, t_]. Cauchy's theorem now yields 

{g,\<f>f) = <P{z{t),Q')f{z{t)). (14) 
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On the other hand, if there is a pole at t = tp away from the cut in the 
complex t-plane (i.e. for t > one would obtain 

Res{(f)f,z{tp)} 



{gtm) = <j>{z{t),Q')f{z{t)) + 



z{tp) - z{t) 



(15) 



The residue can either be approximated or eliminated completely by the 
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bp{z,Q^)=(P{z,Q' 



z - z{tp) 



(16) 



1 - z{tp)z 

where tp is assumed to lie in the range so that 0p is positive for 

z = z{t) with t in [0,t_]. This replacement cancels the pole of f{t) so that 
( [T^ ) holds. The crucial property of this transformation is the fact that, since 
\{z - z{tp))/{l - z(tp)z)\ = 1 ioT z on the unit circle, {(j)f\(l)f) = {4>pf\4>pf) 
and the QCD constraints are left unchanged. 
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Figure 1: Bounds on the form factors f^{t) (on the left side, t increasing to the 
left) and f~^it) (on the right side), derived without any additional constraints 
and plotted over the whole kinematical range of [0, t_] {t in GeV^). 



Because of the positivity of the inner product we have 

V {gMf) {9t\9t) J 



(17) 
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which, by ehminating {4>f\(f)f) with (|TT]), gives us the following bounds on 
the form factors 



\fit)\'<JiQ' 



l-z^it)\<l>izit),Q^)\^ 



(18) 



where we have made use of the fact that z{t) is real for t in [0, 

These Bounds, which are derived using only perturbative QCD and which 
do not depend on any model assumptions or input parameters (they are 
therefore quite loose), are plotted in figure ^ They were calculated using 
Q2 = in J(g2) at order 0(a,) 



If we know the value of the form factor fi 
we can define a matrix Ai 



M 



fm<pf) mgt) mgt.) 

{9t\(t>f) {9t\9t) {9t\9h) 
{9tMf) {9t,\9t) {9h\9t^) 



f{ti) at n points i 



{<Pf\9t.)\ 

{9t\9t„) 

{9tA9tJ 



(19) 



V(^tnl0/) {9tj9t) {9tj9ti) ■■■ {9tj9t„) J 
and obtain, using again the positivity of the inner product 

det > 0. 



(20) 



Since {(f)f\(f)f) is eliminated using (|Tl]) and since all components except f{t) 
of this matrix are known — they are either constants or functions oft — the 
inequality (^) leads to bounds on the form factors f{t). The inclusion of 
some known values (which i.e. are predicted by a model or which one can get 
by lattice calculations, see [M) therefore will give us more stringent bounds 



FUt,Uj,)<fit)<F^j,it,Uji) 

with Flo and Fup calculable functions of t with the parameters ti, fi. 
The upper and lower bounds can be written as 



(21) 



-Pit) ± ^c{Q^) ■ A(t) 
a-0(t,Q2) 



(22) 



The shape of Fup^ioit, Q"^) as functions of t depend on the values for ti and 
In (p2|) , a is a scaling constant, [3{t) gives roughly the shape of the form factor 
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while the squareroot in (^) distinguishes the upper from the lower bound. 
In the squareroot, c(Q^) is a constant in t which depends only on the energy 
scale Q^, and A{t) is a function with zeros in every ti. The reason why the 
function A{t) should behave like this is quite clear: since we fixed the value 
of the form factor at certain points ti, the upper and lower bound should 
coincide in these points (for an exact definition of the functions appearing in 
dH) see 0]). 



3 Improving the Bounds 

The method of including 'fixed points' (i.e. to fix the value of the form factor 
at certain kinematical points) which we summarized in ^ is well known and 
used throughout the literature. However, this formalism can be applied in a 
slightly different way: to include the slope or even higher derivatives of the 
form factor. This is desirable because it is possible to obtain the slope e.g. 
for low t from QCD sum rules or for t near the kinematical end point from 
the chiral limit. 

Consider two fixed points ti and t2 = ti + e with e arbitrary but small. 
We get the coresponding values of /(tj) via a Taylor expansion: 



fih) = /i (23) 

(24) 



/(t2) =/(ti) + e|/(t) 



t=ti 



We can now calculate the bounds depending on these values for ti,t2, fi and 
/2 and take the limit e ^ 0. It is useful to define a function ip{z, Q"^) = 
f{z)(j){z, Q"^), in terms of which we obtain 

_ (1 - z!)iM^ - 2^^i + 4) + - - ^!)) 

± ^ ^'^ , ' X 25 

{l-zzifVT^ 

- (1 + zf) + - zl) - m^{l - zlf 

where zi = z{ti), ipi = 'ip{zi,Q^) and 
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It is obvious how to extend this to higher derivatives (including e.g. the 
curvature of the form factor), but the corresponding equations become te- 
dious. In parts of the numerical analysis presented below the curvature has 
been included. 

It has been observed (see [0, im and in, im) that the inequality 



also yields restrictions for the derivatives of the form factor in the end point 
t_. We shall use this input to restrict the possible values of the form factor 
at t- in terms of the slope and the curvature. In order to do this we rewrite 
([Tl| ) into the form 

^I'^dQ me^^)\'<J{Q'). (26) 

In the case of it is now possible to perform a Taylor expansion around 
Zq = {t = t_) which is convergent in the full unit disc, since /° does not 
have poles or cuts in this region. This is due to the fact that no physical 
intermediate states can contribute. We therefore get the inequality 

CO ^ -| s 2 

j(g')>^'(o) + Eh) ^^"^'(0)- (27) 



n=l 



where ■?/'*•"''* (0) denotes the nth derivative of ip at the point z = 0. Note that 
the form factor, z{t) and thus also ip are real in the region < t < t-. 

The form factor is not analytic inside the unit disc, since there is a con- 
tribution of the B* in this channel. Hence one expects that the Taylor expan- 
sion for this form factor converges only within the circle < \z\ < |2;(m^*)|. 
The integral in (|26| ) runs over the unit circle and thus the integration con- 
tour lies outside the radius of convergence of the Taylor series. In order to 
take into account the i?*-pole one has to expand in a Laurent series or to 
subtract the pole. However, closer inspection reveal that these two methods 
are equivalent. We choose to subtract the pole and thus define a function 

m = m - MZiM^^ (28) 

Z Zrj 



and obtain 



2ti Jo 271 e*® — e*®p 

+r'^w)? (30) 



Res^{0/,zp} r^- <^Q_ 

I- Zl Jo 271 



p 
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where we define ip = f(j). Note that the terms hnear in ifj and l/(e*® — e*®'') 
vanish, since only terms without G dependence contribute to the integral. 
Thus we obtain instead of ( |27D 

jm - ^^^MiM > ^2(0) + £ (l)'v5(")^(o). (31) 

i \n./ 




Figure 2: Possible combinations of point, slope and curvature at the end 
point f'^{t^), allowed parameter triplets lie inside the ellipsoid. To clarify 
the position of the ellipsoid, it's contours were drawn on three faces of the 
cube. 

If one redefines J{Q^) to include the residue and if one uses the analytic 



funtion /, (|3TD is the same as (p7|) 



Since (^) is a sum of positive terms, we can cut off the sum at some 
value of n and get (e.g. for n < 2) 

mm? + m)W + Imrwr < aq') m 
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where we re-substituted ip. 

Such equations give us elhpsoids in the parameter space of value, slope 
and higher derivatives of the form factor at the kinematical end point, where 
all allowed parameter combinations lie inside the ellipsoid. If we take, for 
example, < 2 as in (^) we would get a three-dimensional ellipsoid in 
point-slope-curvature-space (see figure 

The resulting constraints are not as good as the ones that can be obtained 
for B ^ D transitions [0, |15[; this is due to heavy quark symmetries, which 
are not as useful in i? — >■ vr decays. 

However, one may exploit the chiral symmetry for the light degrees of 
freedom as an additional physics input. These symmetries hold at small 
momenta of the pions and are therefore applicable in the end point region 
where the momentum transfer to the leptons becomes maximal. One may 
combine heavy quark and chiral symmetry as in |T^-[0 and compute the 
form factors in this limit. One finds [p^Bl 



fit) 



2/. 



1 + 



1 + 



M^ + m'^ + 2AM - 1 



m 



t) 



M'^ + m? + 2AM - 1 
t 



(33) 



+ 



M2 -m2 



1 - 



g{?>M^ -m^ + t) 
M"^ + w? + 2AM - t 



which results in the residue 

Res{/,tp} 

at 



^^gM{M + A) 

Jtt 



t - M"^ + + 2AM 



(34) 



(35) 



where A = m^. — nis ■ The additional parameters which appear in these re- 
lations are the chiral coupling constant g which describes the 55*7r-coupling 
and the ratio of the decay constants fs/ fn- 

The form factors ( P^ ) are valid in the chiral limit, which holds only in 
a small piece of the kinmatically allowed region. In order to extend this to 
the full range, a variety of ansatze have been invented to extend this pole- 
behaviour to small t 



[see e.g 
ov Vub model dependent. 



p!9|- p2|); however, this makes the extraction 
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We will now use the improved unitarity bounds and the chiral limit to 
derive bounds on the form factors. These bounds have the advantage to be 
model independent, however, as we shall see, they are not very tight over the 
whole range of t. This is mainly due to the fact that the parameters g and 
Ib/ fn suffer from large theoretical uncertainties. 




Figure 3: Possible combinations of point, slope and curvature at the end 
point using the chiral limit. 

In the following we shall consider only; here we can avoid the prob- 
lem of the large uncertainties in g and fs/ f-K by considering fractions like 
f+'{t^)/f+{t^) and f+"{t_)/f+{t_) in which the /s //^-dependence drops 
out and the remaining (^-dependence is small, since terms involving g are 
suppressed: 

r'(t-) M+A 1 

aT^ ' (36) 



2M(M-m)(A + m) 1 



g(M-m) 
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= I I (37) 

/+((_) 2M!(Af-m)(A + m)l + ji^- ^ ' 

Thus in the chiral hmit the ratios f^'{t-)/ f^{t-) and f^"{tS)/ f^{t-) are 
practically fixed and we can express the slope s and the curvature c in terms 
of the value p of at t_. This would yield a straight line in a p, s, c plot; 
however, varying the coupling g in a, generous range between 0.25 and 0.5 



23| and allowing for derivations of the chiral limit of about 15%, we obtain 
the rectangular cone shown in figure |^. 

In case of the residue it is not possible to get rid of the g- and fsj J-r- 
parameter, so one has to take the full uncertainties into account when calcu- 
lating the bounds. One could also choose to drop the residue completely and 
work with the redefined (^p{z, Q'^) of (|16D but we decided to use (|3^) because 
the resulting bounds are more stringent. 

Using (|33D and their first derivatives at the end point as input parameters 
in (^) (the extension to higher derivatives is trivial) we get upper and lower 
bounds for < t < t_ which depend strongly on the chiral parameters g 
and fs/ f-K- In addition to that, a combination of figure |^ and ^ can be 
used to constrain the allowed range for the end point value of /(t-) and its 
derivatives. 

The combination of unitarity bounds, extended to include derivatives, 
and the chiral limit now allows us to calculate bounds for the form factors 
/+(t) and fit). 

4 Numerical Results and Discussion 

Inserting ti = t_ {zi = 0) in ( P^ one gets 

^«p,/o = (^1 + ^[z) ± -^=i=\/j-^2_[^,p (38) 



z 



which leads to the bounds 

<p{z,Q^) 0(^,Q2)V i-z^ • ^ ^ 

We used = for the QCD calculations and inserted all possible combina- 
tions of end point value and slope into (p9). We also allowed variation of g 
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and fe/fn from 0.25 < ^ < 0.5 and 1 < fe/U < 1-7 O Within this 
variation we determined the bounds for fixed parameter values and finally 
took the loosest ones as our result. 

We scanned the parameter space of possible ps-pairs by dividing the al- 
lowed p-range into small, equally spaced intervals. For each of these possible 
p-values, there exists an allowed s-range which is also divided into small 
intervalls. A similar procedure was applied for the chiral input parameters. 



12.5 




-2.5 



Figure 4: Bounds on the form factor f'^{t), derived by using a fixed end point 
and varying the slope according to the chiral limit 

Figure ^ shows the resulting bounds for the form factor /"*". The values 
of the form factor at t- vary between 3.2 and 13.5 due to the uncertainties 
in fB/f-rr and g. The two solid lines represent the result for the maximal 
value of the form factor at t_, where the upper (lower) line corresponds to 
the minimal (maximal) possible value for the slope at t_. Similarly, we have 
plotted the corresponding lines for intermediate values of the form factor at 
t- down to the minimal value, represented by the long-dashed line. 
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The upper and lower bounds coincide at t = due to the kinematical 
constraint (0). We have computed /° in the same way as we did for /+, 
scanning the parameter space for p and s at t = t_. The resuhing bounds for 
/° at t = are much tighter that those for f~^{0) and hence the kinematical 
constraint has a significant effect on /+. In order to be conservative we have 
varied the parameters for both form factors independently; that is we did 
not use the correlation between the form factors implied by the chiral limit, 
namely that at t = t_ they are both given by the same two parameters Jb / fir 
and g. 

In figure ^ we have combined the results for and /°, plotting the 
upper bound for the maximal values of f~^{t-) and f^(t_) with the minmal 
slopes and the lower bounds for the minmal values of f^{tJ) and f^{tJ) 
with the maximal slopes. Compared to the standard method (see figure |lD 
the inclusion of the slope and the chiral limit has significantly improved the 
bounds. We have to point out, however, that not any curve within these 
two bounds is allowed as a form factor, since QCD implies relations between 
slopes and the form factor values. This can be seen in figure ^. 








25 20 15 10 5 5 10 15 20 25 

t 

Figure 5: Bounds on the form factors P{t) (on the left side, t increasing 
leftward) and f^{t) (on the right side, t increasing rightward) derived with 
restriction on possible ps-pairs at /(t-) and including the kinematical con- 
straint, plotted over the whole kinematical range of [0,t_] {t in GeV^). 
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Finally one can also use model input into the machinery of unitarity 
bounds. In particular, a model can be used to obtain points and curvatures 
away from t_. 

In figure |] and |^ we have used the ISGW II model 



2^ as an input. 



Fig.p shows the bounds obtained by the standard formalism using two points 
(t = 10 GeV^ and t_) from this model, while in fig.|^ we use these two 
points but also the derivative at t_, obtained again from ISGW II, as input 
parameters. The dashed line is the model prediction itself. 

Another commonly used model is the BSW model . In fig.|] we use 
this model to determine the value of the form factor /"*" at t-, while in fig.^ 
also the slope at t_ was taken from the model. It is interesting to note that 
if one also includes the curvature of /"*" at t_ in the BSW model — as in 
fig.p!0| — the lower bound shifts significantly upwards, such that the model 
becomes inconsistent with the bounds. This indicates that the curvature of 
the BSW model is too large in the end point, since the lower bound comes 
out to be quite high; its value at t = is about one, which is in contradiction 
not only with the BSW model, but also with sum rule calculations [I 
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Figure 6: Bounds on /° (left) and (right) derived by using the ISGW II 
model and using the values at t = 10 GeV^ and t = t_ as input parameters. 
The bounds are drawn solid, the dashed line corresponds to the model. 
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Figure 7: Same as in figure | but including the slope at t = t_ as well. 




5 10 15 20 25 

t 



Figure 8: Bounds on /+ derived by using the BSW model and using the 
value at t = t_ as input parameter. The bounds are drawn solid, the dashed 
line corresponds to the model. 
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The pole form of the BSW model is motivated by the form factor in the 
chiral limit and thus it is no surprise that they have similar curvatures of 
at t_. We take this as an indication that one should not exploit the chiral 
limit to obtain derivatives of the form factors beyond the slope; going beyond 
the slope requires to take into account higher orders in chiral perturabtion 
theory. 

5 Conclusions 

Using unitarity and analyticity to derive bounds on form factors using input 
from perturbative QCD has become a widely used tool, since this method 
allows to constrain form factors over the whole range of in a model inde- 
pendent way. 

This method has attracted some attention in B decays and has been 
combined with heavy quark symmetry; for the b ^ c transitions this has lead 
to stringent constraints on the form factors, i.e. the Isgur-Wise function. 

In B ^ 71 transitions heavy quark symmetries are not as useful and it is 
necessary to obtain model independent information from other sources. The 
general bounds in B —>■ it decays are not very tight and hence also not very 
useful. Including points where the form factor is known from other sources 
improves the bounds; this method has been combined with lattice data to 
obtain model independent bounds for the form factors for i? — > vr transitions. 

In the present paper we have focussed also at the B ^ 71 decays and 
improved the bounds on the relevant form factors one more step by including 
slopes and higher derivatives of form factors, which in some cases may be 
known from other sources. Again from unitarity the value, the slope, and 
even higher derivatives of the form factor at t_ are constrained to lie inside an 
ellipsoid which can be used as an input into the machinery we have proposed 
here. 

Still this does not tighten the bounds very much and some more physics 
input beyond perturbative QCD is needed. Close to the kinematical end 
point t_ chiral perturbation theory is valid and we used the chiral limit of 
the form factors as an input. With this we arrived at bounds which are much 
tighter than the general bounds obtained without any knowledge on the form 
factors. The model independent results of our paper are shown in fig.§ and 
fig-l- 

Finally one may also combine models with the unitarity bounds from 
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QCD. Taking points, slopes and curvatures from these models one may test 
the consistency of a given model with QCD. An extensive analysis of the var- 
ious models is beyond the scope of the present paper; we only considered the 
ISGW II and the BSW model as examples. Although both models lie within 
the bounds given in fig.^] (the ISGW II touching the lower bound at this 
still does not mean that they are consisten with QCD; the curvature of the 
BSW model at t_ turns out to be incompatibel with the QCD constraints. 
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